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Abstract 

Analyzing a capillary minimizing problem for a higher-dimensional extended fluid, we find that 
there exist startling similarities between the black hole-black string system (the Gregory-Laflammc 
instability) and the liquid drop- liquid bridge system (the Raylcigh-Plateau instability), which were 
first suggested by a perturbative approach. In the extended fluid system, we confirm the existence 
of the critical dimension above which the non-uniform bridge (NUB, i.e., Delaunay unduloid) serves 
as the global minimizer of surface area. We also find a variety of phase structures (one or two cusps 
in the volume-area phase diagram) near the critical dimension. Applying a catastrophe theory, we 
predict that in the 9 dimensional (9D) space and below, we have the first order transition from a 
uniform bridge (UB) to a spherical drop (SD), while in the I0D space and above, we expect the 
transition such that UB — * NUB — > SD. This gives an important indication for a transition in the 
black hole-black string system. 



The most fundamental theory of nature, which is now believed to be string/M theory, predicts that 
spacetime dimension is higher than four. In four dimensional (4D) spacetime, we have a uniqueness 
theorem of a stable black hole (BH). In higher dimensions, however, in addition to a BH, we find a 
new phase of black objects such as a black string (BS) and expect a transition between those objects. 
Its property may depend on the dimension. Hence, in order to understand the physics in higher 
dimensions, it is indispensable to reveal the mechanism of the transitions between the black objects. 

In a spacetime in which one spatial dimension is compactified on a circle, for example, there are 
three possible black "objects", i.e., a localized black hole (LBH) [1,2], a uniform black string (UBS), 
and a non-uniform black string (NUBS) [2-5]. Although the transitions between them may occur, our 
knowledge about those dynamics has been so far very limited. The UBSs are unstable if the size of 
the compact dimension is larger than some critical length scale, as shown by Gregory and Laflamme 
(GL) [6]. Its fate is, however, still controversial [7]. In 5D and 6D spacetimes, "complete" phase 
diagrams are available [2,4] (see [8] for reviews). From these phase diagrams, we find that the UBS 
transits suddenly to the LBH as the compactification scale gets larger than a critical size [3,9]. On 
the other hand, for the case of the large spacetime dimension {D > 14), we have another expectation 
such that the phase transition occurs smoothly, i.e., the UBS is deformed to the NUBS continuously 
[5]. Thus, our current interest is to figure out the complete phase structure and to study whether 
there exists a stable NUBS and how the unstable UBS reaches NUBS. 

To understand the dynamical features of such "black objects" , the membrane paradigm of BHs [10] 
as well as the BH thermodynamics [11] have been sometimes discussed. Among them, there was a 
remarkable discovery by Cardoso, Dias and Gualtieri [12] that the GL instability and the Rayleigh- 
Plateau (RP) instability, which had been found in extended fluids [13], share many common features 
such as the dimensional dependence of the critical wavenumber and that of the growth rate. They 
also pointed out that there may exist a critical dimension [12] similar to that of the BH-BS system [5]. 

For the stability and dynamics of the extended fluid system, there have been extensive studies in 
the case of the usual 3D space. The equilibria of the liquid bridge/drop system are obtained by the 
variation of the surface area while keeping the volume fixed, which is called the capillary minimizing 
problem (or Plateau's problem) [14]. The stability and dynamics can be investigated by the surface 
diffusion equation [15], which does not increase the surface area [16]. The spherical drop (SD) is 
stable [17,18], while the non-uniform bridge (NUB), which is called the Delaunay unduloid [19], is 
always unstable [20,21]. The uniform bridge (UB) becomes unstable when the linear dimension is 
longer than its circumference [17] (the RP instability). The dynamical collapse of a long liquid bridge 
results in a (self-similar) drop formation [22,23]. 

Here we generalize the capillary minimizing problem to the higher-dimensional spaces in which one 
dimension is bounded. It is a mimic of higher dimensional space with one compactified dimension. 
By extremalizing the surface area while keeping the volume fixed, we obtain axisymmetric liquid 
equilibria non-perturbatively. Then, we construct the phase diagrams and compare the results with 
those of the BH-BS system. We find that the system of SD-(N)UB and the BH-BS system exhibit 
remarkable similarities. Note that we should minimize the liquid surface area, while maximize the area 
of event horizons, which is proportional to the BH entropy. Climbing up the "ladder" of dimensions, 
we find that there exists a critical dimension above which the NUB can be a global minimizer of the 
area. In addition, we obtain a very interesting non-trivial structure of phase diagram near the critical 
dimension, which may tell us how the phase transition appears. 

We consider an axisymmetric hypersurface (liquid surface) in (n+2)-dimensional flat space (n > 1), 
in which one dimension in the direction of the symmetric axis is bounded as z £ [— L/2, L/2], where z 
is the axis coordinate. The surface is given by r = r(z), where r represents the radius of the n-sphere 
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at z of the hypersurface. Then, the volume and surface area are given by the functionals of r{z) as 
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where £l n = -n^ n+l ^ 2 /T{{n + l)/2 + 1) is the volume of a unit n-sphere, and a prime denotes the 
z derivative. The equilibrium configuration will be obtained by extremalizing the action I[r] := 
A[r] — kV[t], where k is a Lagrange multiplier (which is geometrically the constant mean curvature of 
the hypersurface). Thus, we have the following Euler-Lagrange equation, 
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where the multiplier k is determined by keeping the volume constant. This equation has two trivial 
analytic solutions; UB, i.e., the uniform bridge with the radius r = ro and k = n/ro, and SD of 
r 2 + z 2 = Rq, i.e. (n + 1)-D sphere with the radius Ro and k = (n + 1)/Rq. 

The static perturbations around the UB gives a simple indication for the similarities between the 
present system and the BH-BS system as shown in [12]. By setting r = ro + er\{z) and k = n/ro + eK\ 
with e C 1, and solving the linear-order perturbation equation r" + nri/rQ = — K\ with the boundary 



conditions of J"]Jz=o = and r'-A^i^ = 0, we find n = r| ^ cos(^/nz/r ) , where r^ 1 is an integration 
constant, and K\ = from the volume constraint. This static solution gives the critical mode of the 
RP instability. We find the critical radius of UB as r cr = y/nL/2ir below which the UB is unstable. 
This sinusoidal onset mode and its dimension dependence [12] are quite similar to those of the GL 
instability [24]. 

Now we show the existence of a critical dimension beyond which the NUBs become stable. The 
present system with ([1]) has a scale invariance, that is, we find one parameter family of solutions by 
rescaling such as L — > PL, r — > (3r, z — > (3z. So we introduce a reduced volume and an area by 



H := V/L 
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which are invariant under the rescaling of length. We denote the quantities associated with the critical 
uniform bridge by the subscript Note that a for a spherical drop, osd, keeps constant as 

\i changes. We have numerically obtained one parameter family of NUB solutions, using the first 
integral of Eq. ([2]), i.e. 



1 



r' 2 + U{r) = , U(r) :-- 



1 



(n+1 



)2 r 2n 



« 2 (r n + 1 + C) 2 



(4) 



where C is an integration constant, and introducing the parameter a := r_/r + which is the ratio of 
the smallest radius (r_) of NUB to the maximum one (r+). A series of a may describe a transition 
from the UB to the SB because a = 1 and correspond to the UB and SD, respectively. 



First, we show the results for two typical dimensions; 4D and 12D, which are representatives for lower 
dimensions (1 < n < 6) and higher ones (n > 10), respectively. The transition range of dimensions 
is 7 < n < 9, which we will discuss later. Several equilibrium configurations for 12D embedded in a 



2 



H/Ha = 1.00 Li/(i CI = 0.980 iij Ha = 0.676 n/n CT = 0.309 \ij\i cx = 0.279 




Figure 1: A uniform bridge (UB), non-uniform bridges (NUBs), and a spherical drop (SD) at the transition in 12D. 
(j, and fj, CT are the reduced volume defined in the text and its critical value. r cr is the radius of the critical uniform 
bridge. The transition occurs from the left to the right as L increases while keeping the volume fixed. The SD appears 
at fi/ncr = 0.279 (< fitouch / fJ-cy = 0.330) via the first order transition. 



schematic 3D space are shown in Fig. [TJ We also show the volume-area diagrams for 4D and 12D in 
Fig. [2J In both dimensions, as a decreases to zero, the NUB reaches the SD which "touches itself" 
at the boundary z = ±L/2 when ji = /U toU ch- For the 4D space, we find that fx increases as the NUB 
smoothly approaches the critical drop, and the NUBs exist only in the larger area region compared 
with those of SDs and of UBs. Hence we suspect that these NUBs are unstable like those in 3D case 
[20,21]. This will be justified when we adopt catastrophe theory (see later). For the 12D case, however, 
the NUB appears in the region with smaller area than those of SDs and of UBs, and it approaches a 
critical drop with a small cusp structure (see Fig. [2]). We expect these NUBs (except for the upper 
branch of a cusp structure) are stable because the areas are smallest for a given volume V and length 
L. This is also confirmed by the argument of catastrophe. This difference depends on whether the 
volume-area curves of SDs and UBs intersects (4D) or not (12D). 



Adopting the least area law, we can predict the most favored (and probably stable) state from 
the volume-area diagrams in Fig. [2j Consider the case fixing the volume and stretching the length L, 
i.e., decreasing fx. In 4D, because the area of NUB is always larger than others, as ll decreases, the 
favored state shifts from the UB to the SD at their intersecting point I with a discrete transition of 
the configuration. This phase diagram is qualitatively similar to the mass-entropy diagram of the 6D 




Figure 2: The volume-area diagrams for 4D (left) and 12D (right), including the uniform bridge (UB), non-uniform 
bridge (NUB) and spherical drop (SD). We enlarge the diagram near the SD in 12D, where we find a cusp structure. 
The solution with the smallest area for a given volume is energetically favored. Hence the UB in 4D and the NUB in 
12D transit to the SD at the intersecting points, I and I', respectively. 
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Figure 3: The tension-volume diagrams for 4D (left) and 12D (right). In 4D, the UB transits to the SD with a finite 
jump of A as the arrow points, while in 12D, the transition from the UB to the NUB is continuous. The transition to 
the SD is, however, discontinuous although the jump is very small (see Fig. 



BH-BS system (Fig. 6 in Ref. [2]). Note that there are a few differences, that is, one is the physical 
role of the area mentioned before, and the other is that the SDs do not deform unlike the LBHs, 
which stems from the neglect of gravity. While in 12D, the favored state shifts from the UB to the 
NUB successively, and eventually transit to the SD. This transition to the NUB is of second order 
or higher. However, there exists a small cusp near the SD. So the transition to the SD occurs at the 
point I', where the two curves of NUB and of SD intersect. This transition is of first order. We expect 
that there also exists this type of cusp in the higher dimensions, although we could only confirm it 
numerically until 15D because of difficulty of numerical accuracy. 

Here we introduce an order parameter of the present system in order to compare our phase structure 
with that in the gravity side. We find that the following parameter, 

L ( dA\ . ^uda . . 

A:= aUi), = - ( " +2 ^- (5) 

which we call the reduced tension, can be an appropriate order parameter. This A is analogous to the 
relative tension in the BH-BS system [26]. A is constant for UB as Atjb = l/(n + 1), and it vanishes 
for SDs. For NUBs, A is evaluated from the gradient of the Li-a curve. In Fig. [3j we show (A, li) 
diagrams. We again find a qualitative agreement between the (A, li) diagram for the 4D space and 
the tension- mass diagram of BH-BS system in the 6D spacetime (see Fig. 5 in Ref. [2]), which had 
been originally conjectured in [9]. As the control parameter li decreases, the UB transits to the SD 
at the intersecting point, li = \x\ (> Li cr ), with a finite jump AA = 1/3 (= Atjb)- While in 12D, 
the UB transits to the NUB continuously. The transition from the NUB to the SD at the point I' is 
discontinuous with a very small jump AA ~ 1.1 x 10" 7 (< Atjb)- 



Now, we discuss the behaviors in the transition range of dimensions (9D, 10D, and 11D spaces). In 
Fig. HI we show the (li, a) diagram and the three dimensional (A, li, a) diagram for 10D space. We 
find that the NUB has two cusp structures in the (li, a) plane. The cusp near the SD is the same 
as that found in 12D. A new cusp also appears near the UB. The upper branch (the larger area) 
corresponds to the curve of the NUB which appeared in 4D, while the lower branch (the smaller area) 
corresponds to the curve of the NUB found in 12D. Because we find the cusp structure, we can adopt 
the catastrophe theory by setting a, li, and A to be a potential function, a control parameter, and 
a state variable, respectively [25]. We expect that the lower branch is stable while the upper one is 
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Table 1: The expected transition patterns in terms of space dimensions. The first ("second") order & first order 
means that we have the first ("second") order transition for UB —> NUB, and then the first order transition for NUB 
— > SD. The "second" means that it may be of second order but might be the higher order transition. For much higher 
dimensions than 15, we have not confirmed the existence of a cusp because of numerical difficulty, but we may expect 
that it is of first order because the topology change is associated in this transition. 



space dimension (n + 2) 


3 ~ 8 


9 


10 11 


12 ~ 


transition type 


UB — » SD 


UB — > NUB — » SD 


order of transitions 


first order 


first order & first order 


"second" order & first order 


number of cusps in /i-a diagram 


no cusp 


two cusps 


one cusp 




Figure 4: The volume-area and tension-volume-area diagrams for 10D. Although we find two cusp structures in the 
(/x, a) diagram, the curve of the NUB in the 3D phase space is smooth, and connects the UB and the SD at its both 
ends. 

unstable. This expectation is justified from the fact that the area in the upper branch is always larger 
than that of the UB, which means that the NUBs are unstable against small perturbations just as the 
3D case. Such interpretation is consistent with our expectation about stabilities in 4D and 12D. The 
Landau-Ginzburg argument also can be applied [9]. The similar cusp structures are observed in 9D 
and 11D spaces. In the 9D case, however, such a structure appears always in the larger area region 
than that of UBs. So this has nothing to do with a realistic transition. 

We summarize the results in Table [TJ The results strongly support the existence of the critical 
dimension on the gravity side [5] and fills a "missing link" in the phase diagrams of the BH-BS system, 
unknown except 5D and 6D [2,9]. 

Finally, we should stress that all equilibria obtained here are static solutions of the surface diffusion 
equation [15]. Although we could analyze the stability from the second variation of the area [20], one 
can also examine the stability and the liquid bridge/drop transition process by solving the dynamical 
equation. What is the most interesting is whether or not we confirm the local/global convergence to 
the NUBs for n > 8. Such analysis would be a footstep to understand the dynamics of the black 
strings in higher dimensions. 

We would like to thank O. J.C. Dias for suggesting a possible extension of his work, and B. Kol and 
R. Emparan for fruitful discussions. U.M. is supported by the Golda Meir Fellowship, by The Israel 
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Research Cooperative Program, and by the Waseda University Grants for Special Research Projects 
and for the 21st-century COE Program at Waseda University. 

Note added — After the first submission of a manuscript to arXiv, one of the present authors wrote 
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an article [27] in which the linear and non-linear dynamics of the RP instability were investigated in 
the framework of the surface diffusion. In particular, it was shown that the NUB indeed serves as an 
endpoint of the RP instability around the critical dimension and above. 
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